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ABSTRACT 



Context. Galactic black hole coronae are composed of a hot, magnetized plasma. The spectral energy distribution 
produced in this component of X-ray binaries can be strongly affected by different interactions between locally injected 
relativistic particles and the matter, radiation and magnetic fields in the source. 

Aims. We study the non-thermal processes driven by the injection of relativistic particles into a strongly magnetized 
corona around an accreting black hole. 

Methods. We compute in a self-consistent way the effects of relativistic bremsstrahlung, inverse Compton scattering, 
synchrotron radiation, and the pair- production/annihilation of leptons, as well as hadronic interactions. Our goal is to 
determine the non-thermal broadband radiative output of the corona. The set of coupled kinetic equations for electrons, 
positrons, protons, and photons are solved and the resulting particle distributions are computed self-consistently. The 
spectral energy distributions of transient events in X-ray binaries are calculated, as well as the neutrino production. 
Results. We show that the application to Cygnus X-1 of our model of non-thermal emission from a magnetized corona 
yields a good fit to the observational data. Finally, we show that the accumulated signal produced by neutrino bursts 
in black hole coronae might be detectable for sources within a few kpc on timescales of years. 

Conclusions. Our work leads to predictions for non-thermal photon and neutrino emission generated around accreting 
black holes, that can be tested by the new generation of very high energy gamma-ray and neutrino instruments. 

Key words. X-rays: binaries - radiation mechanisms: non-thermal - gamma-rays: general - neutrinos 



1. Introduction 

Some physical processes that take place near a black hole 
can be inferred from X-ray observations of accreting bina- 
ries. There is strong evidence that the X-ray emission in 
some of these systems is powered by accretion onto a black 
hole. Accretion processes around compact objects are mod- 
eled by hydrodynamic equations of viscous differentially- 
rotating fluids. The standard disk model is the most famous 
model based on this set of equations, and was developed by 
[Shakura & Sunyaev (1973) and Novikov & Thorne (1971). 
The model predicts that the accreting gas forms a geo- 
metrically thin and optically thick disk, producing a quasi- 
blackbody spectrum due to thermal emission. The effective 
temperature of the accreting gas is in the range 10^ — 10^ 
K, depending on the mass of the compact object and the 
accretion rate. 

The observed X-ray spectrum, however, is too hard in 
many cases to have been produced by a standard disk alone, 
since the hard X-ray component corresponds to a temper- 
ature of ~ 10^ K. 

To explain the complete X-ray spectrum of Cygnus X-1 
-which is arguably the most well-studied black hole candi- 
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date in the Galaxy- an extra componen t is usually added 
to the disk, the so-called corona (e .g., iDove et all 119971 : 
iGierlinski et al.lll997t lPoutanenlll998D . In this context, the 
corona is coupled to the disk by the magnetic field, and the 
plasma is heat ed by reconnecting magnetic loops emerging 
from the disk (jGaleev et al.lll979[ ). 

As a result, the soft photons emitted by the disk gain 
energy by successive Compton upscattcrings in the corona. 
This process is known as Comptonization, and it is the 
most accepted mechanism to explain the broad-band hard 
X-ray spectra of Galactic black holes. Another feature that 
can be explained by the presence of the corona, is a hard- 
ening of the spectra at ~ 10 keV, which is attributed to 
the Compton reflection of hard radiation (e mitted by the 
corona) from a cold materi al, i.e. the disk (jWhite et al.l 
119881: iGeorge fc Fabianl[T99l . 

Since their ear l y applications to Cygnus X-1 
(jShapiro et al.l 119761: llchimarul 119771 ). different models 
have attempted to explain the complete X-ray spectrum 
of black holes in binary systems. The observational data 
of the spectra are used to set some constraints on the 
geometry of the source. For example, the covering fraction 
of the hot cloud as viewed from the soft photon source 
can be estimated from the observed spectral slopes, and 
the reflection bump limits t he solid angle subten ded by 
the disk around the corona (jPoutanen et al.lll997D . There 
are two specific geometries that seem to reproduce well 
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most of the observed spectra: the disk + corona 'sombrero' 
model and the advection dominated accretion Row model 
(ADAF). 

In the 'sombrero' model, the corona is added to the disk. 
It is usually assumed that the corona can be represented 
by a homogeneous spherical cloud of radius Rc around the 
compact object; the cold disk (with an inner radius nn) is 
truncated at a certain distance from th e compact object 
in the hard stat e of Galactic black holes (jPove et al.lll997t 
lPoutanenlll998D . The ratio r-m/Rc depends on the amount 
of Compton reflection observed, and can be obtained from 
the energy balanc e and the electron /posi tron pair balance 
pove et alJ Il997t iPoutanen et~all I1997D . In the case of 
Cygnus X-1, this ratio is estimated to be rin/Rc — 0.8 — 0.9, 
which means t hat only a short fracti on of the disk is within 
the hot cloud (jPoutanen et al.lll997t ). 

On the other hand, the ADAF model is a self-consistent 
solution of the hydr odynamic equat i ons of viscous rotat- 
ing flo w s (see, e.g.. Ilchimarul 19771 : iNaravan fc Yil 11994 
Il995allbl : lAbramowicz et al.lll995ir The main characteris- 
tic of this solution is that most of the energy is accreted 
onto the compact object and the gas is unable to cool effi- 
ciently, mainly because of the low density. Since the plasma 
is a poor radiator, the viscous energy is stored in the gas 
as thermal energy and the gas temperature becomes high. 
This causes the accreting gas to swell, hence the geometry 
of ADAFs is quasi-spherical. This geometry is similar to 
the geometry of the corona-|-disk 'sombrero' model already 
mentioned above. 

Besides the X-ray spectra observed in X-ray bi- 
naries (XRBs), sources such as Cygnus X-1 produce 
steady emission up to a few MeV (jMcConnell et al.ll2000t 
ICadoUe Bel et aLll2006f) , that is indicative of a non-thermal 
contribu t ion to the spectral energy distribution (SED). 
iLi et al.l (|1996f) treated simultaneously the transport and 
acceleration mechanisms for leptons obtaining good fits to 
the s pectra of Cygnus X -1 observed by COMPTEL. In addi- 
tion, jLT^^iilegllH^) demonstrated that under the phys- 
ical conditions expected in the surroundings of accreting 
black holes, it is possible to accelerate electrons out of a 
Maxwellian distribution, resulting in a non-thermal tail. 

More recently, several works have been devoted to 
studying the effects of non-thermal populati ons of elec- 
trons in Galactic black hole coronae (e.g.. iBehnont et aP 
120081: iMalzac fc BelmontI l2009t IVurm fc PoutanenI 12009^ 
as well as in accre t ing supermassiye bl ack holes (e.g., 
iBelmont et aD l2008t IVeledina etaD l201lD . These works, 
however, have not studied the effects of a hadronic com- 
ponen t in the source . 

InE omero et all (I20inch . it was shown that the pres- 
ence of hadrons can also explain the non-thermal emission 
detected in Galactic black holes. In the present work, we 
show the results of a theoretical study of the effects of the 
injection of non-thermal particles, both electrons and pro- 
tons, in a magnetized corona around a black hole, that has 
many refinements compared to previous works. In particu- 
lar, a self-consistent treatment of photon and particle trans- 
port is now presented. We solve the set of coupled kinetic 
equations for all types of particles, including photons, hence 
the treatment of absorption is straightforward. We attempt 
to estimate the SED produced by these relativistic parti- 
cles, and explain the origin of the non-th e rmal tail observed 
in some XRBs (iMcConnell et al.l I2OOOI: iCadoUe Bel eTall 
120061: iJourdain et al\\201^ . 
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Fig. 1. Schematic geometry of the system in the low- hard 
state. 



The completion of IceCube opens new possibilities for 
neutrino detection, so we also study the neutrino produc- 
tion in the source. 

The structure of this paper is as follows. In the next 
section, we present the basic model used to obtain the val- 
ues of the relevant parameters, and fix the characteristics of 
the medium where the non-thermal particles are injected. 
In Sec. [21 we discuss the treatment of the physical pro- 
cesses, along with the details of the numerical method used 
to solve the corresponding equations. We first study the 
system in a steady st ate, to test our code and reproduce 
the results obtained in iRomero "eTall (l20To3) . In SecH we 
then consider the case of a transient event, including the 
time-dependence of the transport equations. Finally, we es- 
timate the electromagnetic emission as well as the neutrino 
fiux. 



2. Basic model 

In accordance with estimates for Cygn us X-1, we assum e 
a black hole of mass Mbh = 14.8M0 (jOrosz et al.ll201l[ ). 
We study the system in the low-hard state, which is the 
state where the X-ray spectrum is dominated by the coronal 
emission. 

The size of the region where the hard radiation is pro- 
duced, is limited by the variability observed in the spectra 
of Galactic black holes. The minimum variability timescale 
is on the order of milliseconds, hence in the har d state the 
corona lies within '^ 20 — 50 rg ()Poutanenlll998f ). where rg 
is the gravitational radius (rg = GM/c^). We considered 
a spherical corona with a size of R^ = 35rg (see Fig. [1]), 
and assumed that the lumin osity of t he corona is 1 % of 
the Eddington luminosity (Es in et all [1997). which results 
in Lc = 1.9 X 10^^ erg s~^. 

In ADAF models, ions and electrons interact only 
through Coulomb collisions and there is no coupling be- 
tween the two species. In this case, the plasma has two 
temperatures, with the ion temperature (T^ = 10^^ K) be- 
ing muc h higher than the elec tron temperature (Tg = 10^ 
K) fe.g.. INaravan fc Yilll995a l[bl). 

We initially consider that the system is in steady 
state, thus the equipartition of energy between the diffcr- 
ent components o f the system is a reasonable assumption 
(lEsin et al.lFT997l) . The power of the jet observed in the 
low-hard state of X-ray binaries is related to the magnetic 
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field, which in turn is coupled to the corona. We then as- 
sume equipartition between the magnetic energy density 
and the bolo metric photon density of the corona Lc (as 
done by, e.g., iBednarek fc Giovannellil [20011 ) . At the same 
time, since the magnetic field launches the plasma into a 
jet, we also assume equipartition b etween the magn etic and 
kinetic energy densities (see e.g., IZdziars kil ll998D . These 
conditions allow us to obtain the values of the main pa- 
rameters. Specifically, the conditions are 
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where B is the value of the mean magnetic field in the 
corona and Ue, rii are the electron and ion plasma densities, 
respectively. 

The X-ray emission of the corona is characterized by a 
power law in photon ene rgy e with an ex ponential cut-off 
at high energies (e.g. .Romero et al.ll20d2h 

nph(e) = Ve-"e-^/^"erg-i ^^-3^ (3) 

We adopt a = 1.6 and Cc = 150 ke V, which are typical 
of Cygnus X-1 (jPoutanen et al.lll997D . The photon field of 
the accretion disk is modeled as a blackbody of temperature 
fcTd = 0.1 keV. Both the X-ray emission of the corona and 
the radiation field of the disk are considered as seed photon 
sources for Compton scattering and photomeson produc- 
tion in relativistic particle interactions. 

In ADAF models, particles are advected onto the com- 
pact object at a mean velocity of w ~ 0.1c. The accre- 
tion timescale is about a second, which i s considerably 
shorte r than the cooling rate for hadrons (jRomero et al.l 
l2010d ). As a result, only high energy protons are able 
to lose a significant amount of energy, and most of the 
power injected in terms of relativistic hadrons falls into the 
black hole. Since most of the electron/positron pairs are 
produced by hadronic interactions, in advective-dominated 
coro nae the h igh e nergy emission is likely produced in the 
iet (iMarkoff et all I2OOII: IVila fc Romero! 120101: IVila et al.l 

l2om 

In iRomero et all (l2010cD . a static corona model was 
shown to be capable of producing the non-thermal emis- 
sion of Cygnus X-1. In this work, we consider only this 
possibility. In a static corona, which is supported by mag- 
netic pressur e, the relativistic p a rticles can escap e mainly 
by diffusion (jBeloborodovl 119991: lNavakshin|[l999D . In the 
Bohm regime, the diffusion coefficient is D{E) = rgc/3, 
where Tg — E/{eB) is the gyroradius of the particle. The 
diffusion rate is 
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We consider the injection of the non-thermal particle 
distributions of both electrons and protons in this corona. 
The mechanism of particle acceleration in black hole coro- 
nae is likely related to magnetic reconnection, which is es- 
sentially a topological reconfiguration of the magnetic field 
caused by a change in the connectivity of the field lines 
Several works have been published on particle ac celera- 



tion through magnetic reconnection (e.g. Schoppe r et al 



ll998HZenitani k Hoshinoll200lHde Gouveia Dal Pino et al 



l2010t ) The basic idea is that a first-order Fermi mech- 
anism takes place within the reconnection zone caused 
by two converging magnetic fluxes of opposite polar- 
i ty that move toward each other w ith a velocity v^cc 
(|de Gouveia dal Pino fc Lazarian|[2005[ ). The resulting in- 
jection function of relativistic particles is a power-law with 
an index F ~ 2.2. 

The detailed analysis of lDrurvl (|2012l) supports the idea 
that magnetic reconnection can lead to type I Fermi accel- 
eration in a similar way as diffusive acceleration is driven 
by shocks. According to Drury, the spectral index does not 
have a universal value of 2.2; it instead lies somewhere in 
the range 1 < F < 3. Despite this caveat, and given the un- 
certainties in the acceleration models, we adopt the value 
of F = 2.2, which is consistent with all simulations imple- 
mented so far. 

In our model, the acceleration mechanism is not in- 
cluded as a term of the energy gain in the transport equa- 
tion, but used to fix the injection function of primary elec- 
trons and protons and determine the maximum energies 
that relativistic particles can achieve. 

As for standard first order Fermi acceleration, the ac- 
celeration rate t~J^ = E~^dE/dt for a particle of energy E 
in a magnetic field B is given by 



t~ = 

acc 



•qecB 



(5) 



where is a parameter that characterizes the efficiency of 
the mechanism in the magnetized plasma. It is given by 
(lDrurvlll983l: IVila fc AharOTLia3l2009l: Idel Valle et aII[20Tl 
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The recon nection speed in violent re c onnection events i s 
^'rcc ^ VA (jLazarian fc Vishniadll999t iKowal et al.|[20TTl ). 
where va is the Alfven velocity, given by 



52 



VA = 



ATTiripn 



(7) 



In the corona, the Alfven speed is 0.5c, yielding an ac- 
celeration efficiency of 77 ~ 10"^. 

The power available in the ma gnetized system can be 
estimated by (|del Valle et al.ll201l[ ) 



B^ 
Stt 



AVA, 



(8) 



where A ~ AttR^. This yields a power available for non- 
thermal processes of 15 % Lc. The total power injected 
into relativistic protons and electrons, Lrci, is assumed to 
be a fraction of the luminosity of the corona, Lrci = ft oiLc, 
with Qroi < 0.15. The way in which energy is divided be- 
tween hadrons and leptons is unknown. To deal with this 
uncertainty, it is useful to define the parameter a as the 
ratio of the power injected in protons t o the on e injecte d in 
electrons, a = Lp/Le- Following Romer o" et al.l (j2010d ). we 
consider models with a = 100 (proton-dominated case, as 
for Galactic cosmic rays) and a — 1, that is, models with 
the same power injected in both hadrons and leptons. In 
our model, the injection function is assumed to be both 
homogeneous and isotropic. 

We model the escape of radiat i on fro m the region using 
the treatment described in lCoppj (|l992f ) 
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Table 1. Main parameters of the model. 



Parameter 


Value 


Mbu'- black hole mass [AIq\ 


14.8" 


Re', corona radius [rg] 


35 


rin/Rc'- inner disk/corona ratio 


0.9*" 


le- electron temperature [KJ 


10 


Til ion temperature [K] 


10^2 


tc- X-ray spectrum cut-off [keV] 


150 


a: X-ray spectrum power-law index 


1.6 


77: acceleration efficiency 


10"^ 


Be: magnetic field [G] 


5.7 X 10^ 


ni,ne: plasma density [cm~^] 


6.2 X 10" 


kT: disk characteristic temperature [keV] 


0.1 



Notes. Value for Cygnus X-1 l|Orosz et alJl201ll ). 
"■'^ This is the typical value in models where the main source 
of seed photons for Comptonization are provided by the disk 
ijPoutanenlliggi : iHaardt fc Maraschilll993l '). In our model, there 
are other relevant sources of photons that can relax the con- 
dition Tin/Rc = 0.9; however, we only use this parameter to 
constrain the size of the corona, which for Cygnus X-1 ca nnot 
differ significantly from the one adopted in our work (jPoutanenl 
[1991 ). 



Rc 



where 



09 



for X < 0.1, 
for 0.1 < X < 1, 
for X > 1, 



2Rc < (TKN-E'ei > N^±, 



(9) 



(10) 



(11) 



X — E-y/nieC^, ctkn is the Klein- Nishina cross-section, 
N^± is the number density of pairs defined by N^± — 
J dE^±N^±(E^±), and < aKNE^± > represents an aver- 
age over the particle distribution such that photons with 
energy > rrieC^ escape out of the source in a time Rc/c. 
Table [1] summarizes the values of the relevant parameters 
in our model. 



3. Treatment of radiative processes 

Computing the SEDs of black hole coronae, as well as 
other magnetized plasmas, is a complex task, since it must 
include, among other issues, a detailed knowledge of the 
plasma characteristics and the microphysical processes. The 
first method use d to treat this problem was a Monte Carlo 
simulation (e.g., Aharonian et al.l [19851 : IStern etallllQQSt 
iPilla fc Shahamlll997D . The main problem for this approach 
is the small number of high-energy photons, which leads 
to low quality photon statistics. On the other hand, it is 
usually easy to mo del the radiative transfer processes (see 
iPellizza et al.ll20Tol for a three-dimensional code). 

A second method for estimating the spectra of 
compact sources involves solving the kinetic equa- 
tions ('e^g..'Light man fc Zdziarskilll987t[CoDDi fc BlandfordI 
119901: [Coppi 199^ The different interactions of particles 
with the fields of the source ensure that it is unavoid- 
able to study a wide range of particle energies, hence the 



computation of radiative processes is quite complicated us- 
ing this approach. The main advantage of this method is 
that the transport of photons is solved self-consistently. 
Since, in general, these equations are solved numeri- 
cally, the availability of computational resources over 
the past decade has allowed the significant improvement 
of thi s approach (see, e. g., Aharonian fc Plvasheshnikovj 
20031: iMalzac fc Belmont ^09.; Poutancn & Vur m 20091: 
Vurm fc Poutanenll2009i r 

To complete the treatment initiated in iRomero et al.l 
(l20T0^ . we take the second approach. As mentioned in 
the previous section, we are interested in the study of the 
injection of non-thermal particle distributions of electrons 
and protons in the system. Once protons are injected into 
the corona, they interact with both the photon and mat- 
ter fields, producing pions. In addition, the charged pions 
decay producing muons, so we also take into account the 
presence of these transient particles. 

An accurate description of a hot, magnetized plasma 
such as the corona should also treat the processes of 
pair production and annihilation, hence we include elec- 
tron/positron pairs. 

The main channel for secondary pair production in our 
model is photon-photon annihilation. The most important 
background photon field for pair creation is the thermal X- 
ray radiation of the corona. The e nergy spectrum of pairs 
has been studied, for examp l e, by lAharonian et"aLl (|1983D 
and iBottcher fc Schlickeiseil (|l997r ). Under the conditions 



of erg ~^ 



< E. 
,-1 
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the pair emissivity Q 



: (Ee) (in units 
can be approximated by the expression 



77- 



i 77— 



= (£^e±) 



3 cctt 
32 TOgC^ 



4e2 



de-yduj 



n^(e^) nph(a;) 



■1Tc(«7 fc) 



leie-r - le) 



■In 



47^0; (e^ - 7e) 



2{2e^uj - l)e^ ^ ^ 



1 



(12) 

Here 7e = E^jm^f? is the Lorentz factor of the electron, 
= E-y/nieC^, and uj — e/rrieC^ are the dimensionless pho- 
ton energies. 

Another important source of electron/positron pairs 
when protons are present is the Bethe-Heitler process. To 
est imate this contribution, w e use the treatment described 
bv lChodorowski et al.l ()1992[ ). 



3.1. Radiative losses 

Loss terms in our equations include synchrotron radi- 
ation, inverse Compton (IC) scattering, and relativistic 
bremsstrahlung for electrons and muons. We also consider 
photon production by pair annihilation. For protons, the 
relevant mechanisms are synchrotron radiation, photome- 
son production, and hadronic inelastic collisions. 

A complete discussion of cooling times due to syn- 
chrotron radiation, IC scatte ring, and hadronic interac tions 
can be found, f or exam ple, in I Vila fc AharonianI (j2009f ) and 
iRomero et al.1 (|2010c[ ). The injection of secondary particles, 
such as pions and muons, is also discussed in these works. 
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The number of pairs is modified by both creation and 
annihilation. A useful ap proximation to the annihilation 
rate is given bv lCoppi fc Blandford (,1990i) . 



3 O-Tc(TOeC^)^ 



E± 



(13) 



where N^(Ez^) represents the electron/positron distribu- 
tion (in units of erg~^ cm"^). 

In Fig. [21 we show the cooling rates together with 
the diffusion, decay, and acceleration rates of all particle 
species. The value adopted for i?min is twice the rest-frame 
energy of each type of particle. The maximum energy for 
electrons and protons can be inferred using a balance be- 
tween the acceleration rate, given by Eq. [SJ and the cool- 
ing rate. This yields i^^ax ~ 10 GeV for electrons, and 
-^max ~ 10^ TeV for protons. Particles of such energies sat- 
isfy the Hillas criterion, and can be contained within the 
source. 

The cooling rates for the IC scattering and photomeson 
production shown in the figures are the result of the inter- 
actions of particles with two target photon fields: the X-ray 
emission of the corona and the emission from the disk. 

Cygnus X-1 is a binary system with a massive star 
that produces an intense radiation field. This photon field 
could be considered as an additional target for IC scattering 
since it dominates the bolometric luminosity of the source. 
However, its effect on particle energy losses is negligible. 
This can be shown by a simple analysis: in the Thompson 
regime for IC scattering, the cooling rate {t~^) is propor- 
tional to the energy density of the target photon field, de- 
fined by 



(14) 



Uph = / EphUphdEph. 



This magnitude for both the accretion disk and the corona 
X-ray photon fields is on the order of magnitude ~ 10"'^° 
erg cm~^, whereas for the stellar field at the location of the 
corona it is ~ 10'^ erg cm~'^. The photon field provided by 
the star is then an ineffective seed target for IC processes 
compared with the local X-ray field. 

As can be noted from the figures, no unique mecha- 
nism clearly dominates the energy losses for a given par- 
ticle species. For electrons and muons, there are two rele- 
vant radiative processes, IC scattering and synchrotron ra- 
diation, whereas for protons and charged pions hadronic 
interactions are also important. Since IC scattering and 
photomeson production are processes that depend on the 
radiation field and, at the same time, the photon field is 
affected by all the interactions of particles with the differ- 
ent fields, the particle cooling times and distributions are 
strongly coupled with the transport of photons. This point 
i mplies that it is essenti al to improve the model presented 
in iRomero et al.l ()2010c[ ) , solving the system of coupled ki- 
netic equations for all type of particles. 



3.2. Set of coupled equations 

We determine the relativistic particle and photon distri- 
butions solving the set of coupled transport equations in 



the steady state and assuming spatial homogeneity and 
isotropy. The set of kinetic equations are 

a) Transport of electron/positron pairs and protons: 

^ {hiE)N,{E)) + ^ = g,(i?), (15) 

where i = e+, e^, p. 

b) Transport of charged pions and muons: 



dec 



where i — 7r+, tt^, /i^. 
c) Transport of photons: 



= Q.{E), (16) 



^^!^^ = Q^(i?^) + Qe±^^(iV,±,i?^) ^^^^ 

- Q^~f^e±iNj,Ej). 

Here, Ni{E) represents the steady state of each particle dis- 
tribution (in units of erg"-'^ cm~^), h{E) includes all radia- 
tive losses for a given type of particle, tosc is the timescale 
over which relativistic particles escape from the system, 
'•dec mean decay time for transient particles (pions 

and muons), and Qi{E) is the injection function. 

In Eq. (flT)) . the term Q~f{E^) represents photon injec- 
tion due to several radiative processes 

QliEj) — QsynchriEj) + Qlc{Ej) + QTrO^^^{Ej). (18) 

where Qsynchr(-B^), <3ic(^^7), and QTr0^y-y{Ej) give the con- 
tribution from synchrotron radiation, IC scattering, and 
neutral pion decay to photon injection, respectively. 

The process of pair annihilation is another source of 
photons. The corr esponding annih i lation line emissivity can 
be co mputed as (jSvenssonl Il982t iBottcher fc Schlickeiseil 
Il996h 



Qe±(A^7:^7) 



where 



dE^+dEf.- Re± {Ee- , Ef,+ , Ej) 
xN,+ {E,+ )N,-{E,-), (19) 



3 (TTc(mec2)5 

Rp± = - — _o „^ X 



E\E^ 



(t^m) 



(7^m) 



- £^e+ I + 2TOeCV7r \E-^ -E^-\ + 2mec2/7r 



7CM 



r(^7+ +7- - Ey/meC^^ 



(20) 
(21) 



for Ey > i?g+ , E^- or E^ < E^+ , E^- , or in any other case: 



7CM 



(1+7+7- + (7^- 1)1/2(7^ -1)1/^), (22) 

and 74- — Ep+ / rripC^ and 7_ = E^- /rrieC^. We refer to 
iRomero fc Paredei (|201lD and references therein for for- 
mulae on radiative processes. 
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Fig. 2. Relevant radiative losses in a corona characterized by the parameters of Table [T] 



3.3. Numerical method 

We u se an Adams-Moulton method (see, e.g.. iPress et al.l ^ 
Il992f ) to solve the differential equations in Eqs. ([TS|) - ([T^ . 
This is an implicit multi-step integration method that can 
reach higher orders than other numerical algorithms; we 

use in particular a second order method^ 

Fo llowing the scheme described in IVurm fc PoutanenI 
we define an equally spaced grid on a logarithmic 
scale for the energy of particles 



In Ei = In Smin + i ■ AE, i e [0,i^], 



lnE^ = \nEl,^^ 



(23) 
(24) 



We then obtain a system of linear algebraic equations of 
the form 



(25) 



with the boundary condition that Ni^ = 0, which repre- 



losses, whereas particle injection is included in the vector 



hi{Qi + Q2) 
hi{Qi + Q2) 

h2{Q2+Qz) 



\ 



2 ^m — 1 {^Qi_ 





Qm) 



(26) 



where hj = Ej^i — Ej is the energy step. 

We first solve the transport equations and obtain the 
particle distributions. These are used to estimate to first 
order the non-thermal luminosity. Once we know the non- 
thermal photon injection, we solve Eq. 1171 An important 
property of the photon transport equation dTT]) is its non- 



sents A^(i?max) = 0. The matrix Aij contains the particle linearity. This is because the cross-section of photopair pro- 
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duction {Qyy^^±{N^, E^)) depends explicitly on the pho- Table 2. Orbital and stellar parameters. 
ton distribution. We use the approximation discussed in 
iPoutanen fc Vurm (2009), which consists in taking the pho- 
ton distribution from a previous step, j, to obtain the cur- 
rent injection (step j + 1) of electron/positron pairs. Since 
we firstly consider a steady state, the way to solve the pho- 
ton transport equation is then reduced to a simple iterative 
scheme given by 



Parameter 


Value 


M*: Star mass [A^q] 


20 


Ri,: Star radius [Rq] 


17 


Ti,: Star temperature [K] 


3.0 X 10* 


Porb: Orbital period [days] 


5.6 


a: Semi-major axis [cm] 


2.3 X 10^^ 


i: inclination angle [°] 


27 



Q: 



(27) 



The updated photon distribution is then added to the 
backgroundphoton fields (corona power-law plus emission 
of the disk)[3 to compute the IC scattering and hadronic 
interactions. We calculate the radiative losses and injection 
of particles where appropriate. The transport equations of 
massive particles are then solved, and the new distributions 
are used to compute the luminosity to second order. The 
process is repeated until all particle distributions converge 
to a stationary value. 

3.4. Photon absorption in tlie stellar radiation field 

The binary system Cygnus X-1 is composed of a massive 
star and a comp act object. The m assive star is an 09.7 lab 
star of ~ 20Mq (jOrosz et al.ll201l[ ). The orbit of the system 
is circular, with a per iod of 5.6 days an d an inclination of 
between 25° and 30° (|Orosz et al.H20lH) . 

The star produces an intense radiation field that can ab- 
sorb gamma rays by pair creation within the binary system. 
The photon field of the star is anisotropic, because its in- 
tensity depends on the position of the black hole in its orbit 
(but outside the corona). The gamma ray absorption in X- 
ray binaries with a massive co mpani o n star has been stu d- 
ied, for example , by Herterich ( 1974 ). Carraminana (1992), 



iBednarekl (Il993l. 12000.) Bottcher fc Dermer. (.2005.) .Dubus 
(|2006D . IZdziarski et all (|2009t ). and lRomero et all (|2010a^ 



To estimate the observable spectrum and compare with the 
available data, it is necessary to include an appropriate 
treatment of the absorpti on in the stellar field. W e then use 
the approach described in lRomero et al.l (|2010al ). where the 
case of Cygnus X-1 is considered. 

The star has a radius i?* = 1.5 x 10^^ cm, and we assume 
a blackbody radiation density of temperature = 3 x lO'* 
K. The orbital radius is r^rh = 3.4 x 10^^ cm. Table H 
lists the values of t he companion star and orbit parameters 
(|Orosz et al.ll201lD . 

In Fig. [21 we show the opacity map produced by the 
photon absorption in the stellar field, in the relevant energy 
range and along the complete orbit {(p is the orbital phase 
in units of 27r; we note that (p — — 1 corresponds to the 
compact object at opposition, i.e. superior conjunction). 
As can be seen in the figure, except for close to the inferior 
conjunction (compact object in front of the star, = 0.5), 
the emission is completely suppressed by the stellar photon 
field at energies in the range 10 GeV < E < 120 GeV. 
To better illustrate this effect. Fig. H] shows the modulation 
with the orbital phase of gamma ray emission of the corona 




10.5 11 
Log(E/eV) 



Fig. 3. Map of the absorption produced by the anisotropic 
stellar photon field. 




100 150 200 250 

Orbital Phase [degree] 



300 



As we have mentioned, at the location of the corona the 
stellar photon field can be considered as negligible. 



Fig. 4. Modulation of gamma ray emission in the steady 
state at i? ~ 50 GeV by the anisotropic photon field of the 
companion star. 



in a steady state at a given energy. The stellar photon field 
is almost transparent when the compact object passes in 
front of the star {(f) = 0.5), but almost opaque when the 
compact object passes behind (0 = 0) . These results agree 
with those obtained bv lRomero et al.l (|2010al) . 



3.5. Spectral energy distribution 

Figure [5] shows the main non-thermal contributions to the 
total luminosity. The internal absorption is not included in 
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this plot. The synchrotron radiation of electron/positron 
pairs dominates the spectrum for < 100 MeV. We note 
that given the small size of the corona, the synchrotron ra- 
diation below E < 1 eV is self-ab sorbed. All radio em ission 
of the source comes from the jet (jStirling et al.ll200ll) . 

The relativistic protons injected in the corona collide 
with thermal protons through different channels (Ci and (2 
are the multiplicities) 



P + P ^ P 



+ n + TT+ + Ci7r° + C2(7r+ + vr"), 



p + p-^n + n + 2Tr+ + Citt" + C2(7r+ + tt" 
and with the photon field 

p + j^p + Civr" + C2(7r+ + vr"), 



P + 1 



, + 7r+ +Ci7r" + C2(7r+ +7r-). 



(28) 



(29) 



(30) 



(31) 



(32) 



The main electromagnetic result of these interactions is the 
neutral pion decay 



7 + 7. 



(33) 



This is the most relevant source of photons in the high- 
energy gamma-ray band (see Fig. [SJ. 

One assumption of our model is the equipartition of 
energy between the magnetic energy density and the pho- 
ton energy density of the corona. It is then expected that 
the contribution of synchrotron radiation and IC scatter- 
ing to the total luminosity be comparable. In the analysis 
of Fig. [21 we showed that IC scattering and synchrotron 
radiation are the radiative processes dominating the en- 
ergy losses for electrons at low energies. At higher energies, 
the Klein-Nishina effect becomes important, the IC cross- 
section decreases, and this leaves synchrotron as the main 
mechanism causing electron energy loss. The Klein-Nishina 
effect is also responsible for the diminution of the IC radia- 
tion with respect to the Thompson regime; the gamma-ray 
flux is proportional to the number of interactions, hence 
when the IC cross-section decreases, so does the gamma 
ray emission. This explains why in Fig. [S] synchrotron radia- 
tion dominates the luminosity of the source at low energies. 
Nevertheless, the IC emission of electron/positron pairs is 
comparable to the synchrotron radiation at i? ~ 10^^^ eV. 

In Fig. ini we show the total photon flux produced in 
the corona in two different orbital phases, superior and in- 
ferior conjunction. These positions correspond to the max- 
imum and minimum absorption in the stellar photon field, 
respectively. We also show the SEDs obtained for two dif- 
ferent values of the parameter a. We compare our results 
with ob servations of Cygnus X-1 made by COMPTEL 
(|McConn ell et al. 2000), obtaining good agreement. The 
best-fit model is obtained with = 0.02 for a corona 
dominated by protons (a = 100) and q^ei = 0.03 for a = 1; 
both values are significantly lower that the maximum en- 
erg y available for particle acceleration. 

Jourdain et al.l (j2012l ) reported observations of Cygnus 
X-1 made with the SPI instrument onboard the 
INTEGRAL satellite. Despite no significant emission above 
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Fig. 5. Main non-thermal contribution to the luminosity, 
without considering internal absorption. 



1 MeV being detected by INTEGRAL, the upper limits 
are consistent with the non-thermal tail observed at sev- 
eral MeV by COMPTEL, and in close agreement with our 
model. 

The gap observed in the energy range 10^ < E < 10^ 
keV is produced by the internal absorption in the corona 
and accretion disk fields. To quantify this effect, we show 
in Fig. [7] the opacity as a function of the photon energy at 
different depths inside the corona. Given the high values of 
the opacity, the emission is completely suppressed. This re- 
sult is in accordance with the non-detection of Cygnus X-1 
i n a ste ady state by Fermi. As pointed out in lRomero et al.l 
(I2010d) . ah emission detect ed in this energy range sh ould 
be produced in the jet (e.g.. iBosch-Ramon et al.ir2008[ ). 

The absorption in the stellar field partially suppresses 
the high-energy bump at iJ ~ 10^'^^^^ eV, which makes 
it difficult to detect this source using either the MAGIC 
or VERITAS Cherenkov telescopes. The high-energy emis- 
sion may be detectable by future instruments with higher 
sensitivity and wider energy ranges, such as the Cherenkov 
Telescope Array (CTA). 



4. Transient episode or flare 

It is well-known that X-ray binaries undergo transient ra- 
diative episodes. In 1999, the BATSE instrument detected 
an increase in the lumino sity above 50 keV in Cygnus X-1, 
of an order of magnitude (jStern et al.ll200lD . Moreover, be- 
tween 1995 and 2003, seven outbursts have been reported at 
the location of this source with a signific a nce of 3g or more 
(M azets et al.l 119961: iRomero et alll2002l: iGolenetskii et al.l 
i2003l ). The luminosities above 15 keV of the outbursts were 



in the range 1 — 2 x 10 erg s , which are much higher 
than the typical thermal luminosity in the hard state. 

More recently, the claime d 4.2(7 detec t ion o f Cygnus 
X-1 during a flaring state bv lAlbert et al.l (I2007D and the 
suggested detection by the AGILE satellite ( Tavani et al.l 
2009) constitute the flrst presumed evidence of very high- 
energy gamma ray emission produced around a Galactic 
black hole. 
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Fig. 6. Final flux in a corona + disk characterized by the parameters of Table [TJ We include the 5cr sensitivities for 
different instruments (50 hours of direct exposure for MAGIC and CTA and 1 yr survey mode for Fermi). 
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Fig. 7. Internal absorption due to photopair production in 
the soft photon field of the corona and the accretion disk. 



Another example is the X-ray binary Cygnus X-3, from 
whic h four gamma-ray f lares -were detected by AGILE satel- 



lite (jTavani et al.ll2009(). Variable emission -was indeed de- 



tected by Fermi ( Abdo et 311120091). This emission, ho"wever 
is li kely produced in the jet (|Bednarei3l2010l : lAraudo et al.l 
^0 ; Cerutti ct al. 2011). 

With the aim of studying transient events in the corona, 
"we apply our model to the case of a non-thermal flare, i.e. 
a flare produced by changes in the total power injected into 
relativistic particles (groi)- The dynamical element during 
these type of outbursts is the magnetic field; a sudden injec- 
tion might be the result of both fast and large-scale recon- 
nection events. This assumption is supported by observa- 
tions of solar flares that suggest that magnetic reconnection 
can trigger diffusive accel eration without the requirement of 
strong shock formation (jTsuneta fc Naitdll998l : lLinll2008l: 



IKowal et al.ll201ll) . The overall thermal luminosity can re- 
main constant during these episodes. 

4.1. Particle injection 

Although the light curves of Galactic black holes during 
outbursts or transient episodes can vary from source to 
source, there are some common features. T he rise time 
tends to be much shorter than the decay time (| Grove et al.l 
I1998D , thus the light curves are usually called FRED (Fast 
Rise and Exponential Decay) . A simple analytic expression 
that c an represent this behavior is given by ((Romero et all 
I2010bn 



QiE,t) ^QoE-'^e 



-a~E/E„ 



t/r 



arctan 

2 V Td, 



^(1-e 

'plat 



(34) 



where Tj-isc, Tdoc, and Tpiat are the rise, decay, and plateau 
duration, respectively. We adopt Triso = 30 min, r^ec = 1 
h, and Tpiat = 2 h for a rapid flare. The power-law has the 
standard index of a — 2.2. The normalization constant Qq 
can be obtained from the total power injected into rela- 
tivistic protons and electrons, Lj-oi = Lp + Lg. This power 
is assumed to be a fraction of the luminosity of the corona 
-^rci = QreiLc- As mentioned in the previous section, in the 
steady state the best fit to the observations is obtained 
with 5i-ci = 0.02 for a — 100. During the flare, the number 
of relativistic particles increases. In our model, the power 
injected into the flare doubles that injected in the steady 
state, but larger flares are quite possible, as observed in 
the Sun (Lin 2008). It is assumed that the thermal corona 
remains unaffected during the event. 



4.2. Spectral energy distribution 

In Galactic black hole coronae, cooling timescales are sig- 
nificantly shorter than the flare timescales, which are typ - 
ically of hours or even days (jMalzac fc JourdainI |2000[ ). 
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Fig. 8. Evolution of the luminosity during a moderate non- 
thermal flare. We adopt a — 100. Only the absorption in the 
corona photon field is considered here, since the transport 
of photons is solved self-consistently. The absorption in the 
stellar field, however, is not included in these plots. 

Consequently, the transport equation could be equally 
solved assuming a steady state, and considering changes in 
the luminosity as the flare evolves. Since one of the aims of 
this work had been to develop a code that can be applied to 
different environments, we include the time-dependent term 
in the kinetic equat ions. The transport equations then have 
the following form (jGinzburg fc Svrovatskii|[l964|) 



dN,{E,t) 



dt 

where as before 



_d_ 

dE 



(b{E)N,{E,t) 



N,{E,t) 



h{E)^ 



dE 
'dt 



(35) 
(36) 



We solve the set of coupled equations using the treatment 
described in Sec. 13.31 but now including the time depen- 
dence so the computing time significantly increases. 

Figure \E\ shows the evolution of the electromagnetic 
emission during a day. In this figure, we do not include 
the absorption in the photon field of the star, since it de- 
pends on the orbital phase where the fiare is produced. For 
this purpose, in Fig. |9] we show the absorption coefficient 
related to the stellar field. We estimate the opacity for flares 
occurring at different orbital phases, and we conclude that 
for some values of 4> the absorption of the star is almost 
negligible. For example, flares at energies above 10 GeV 
can be detected at phases </> ~ tt, with instruments such as 
CTA. 



4.3. Neutrino emission 

The electromagnetic flare may be absorbed above E > 10 
MeV. This is not true for a neutrino burst. If a sudden injec- 
tion of relativistic protons occured, the neutrino flux pro- 
duced could be detectable by instruments such as IceCube. 

We consider neutrino production by two main channels: 
charged pion decay 



TT^ ^ fi^ + (37) 

and muon decay 

/i"^ ^ -f F^(z/^) + l^eiVe). (38) 



Thus, the total emissivity of neutrinos is (jLipari et al.ll2007l : 
iRevnoso fc Romerol[2009t ) 



]u{E,t) = Q {E,t) + Q^^,{E,t), 



(39) 



where 



Q^^,{E,t)^ / dE^t-\^^{E^)N^iE^,t)-, 
Je 

^ 6(1 -r^-x) 
E^{l-r^) ' 

with X — E/Et^, = ("i^/m^)^ and 



(40) 



4 r 

Q^^AE,t)^Yl / 



dE. 



p ,-1 



E„ 



ME,)N,AE,,t) (41) 



1=1 ' 
5 „ 9 4 . 



In this latter expression, x ~ E/E^, M{i,2} = /^l 

M{3.4} = Mr '^^ ■ Our calculations take into acco unt pion 
and muon losses, as in IRevnoso fc Romerol (|2Q09t) . 

The differential flux of neutrinos arriving at the Earth 
can be obtained as 



d<P, 
dE 



(42) 



This quantity, weighted by the squared energy, is shown in 
Fig. [TOl which also shows the IceCube sensitivity for one 
year of operation. Assuming that the duty cycle of flares in 
Galactic black holes is around 10 %, the IceCube detector 
will be able to de tect neutrinos from a source at ^ 1.8 kpc 
(iReid et al.ll201lD after ten years of observations. 

Nonetheless, a variability search with carefully binned 
time spans migh t yield positive results long before 
(|Vievro et al.ll2012l in preparation) , avoiding the disadvan- 
tage of the smoothing cause by the averaging of the obser- 
vations. 



5. Discussion 

We have developed a model to deal self-consistently with 
the non-thermal emission from a magnetized corona. The 
assumption of equipartition among the magnetic energy 
density, the bolometric photon energy density of the corona, 
and the kinetic energy density of the plasma is used to esti- 
mate the value of the most relevant parameters. This is sup- 
ported by the X-ray binaries spending a signiflcant amount 
of time in the low-hard state, permitting in turn a signif- 
icant amount of time for field-particle interactions. Under 
this assumption, the model presented in this work is self- 
similar, because if a value of one parameter is changed, the 
other parameters can be re-scaled, producing no significant 
differences to the results and predictions. 
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Fig. 9. Changes in the opacity to gamma-ray propagation for flares produced at different orbital phases. All flares a 
assumed to last ten hours, this is ^ 7.5 % of the orbital period. 
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On the other hand, we found that if we consider a sub- 
equipartition magnetic field, the radiative output can be 
very different. The role of the magnetic field in the cooling 
of the different types of relativistic particles in the corona is 
very important. If lower values of the magnetic field are con- 
sidered, the main radiative losses will be produced by the 
interaction with matter (e.g., relativistic bremsstrahlung) , 
hence the spectra of relativistic particles will be modified. 

The hadronic content of the plasma, given by the 
parameter a, is unknown. We have adopted two values: 
a ~ 100, which represents a proton-dominated corona, and 
a = 1, i.e. a corona with equal contributions from both 
protons and electrons. Both models are capable of repro- 
ducing the non-thermal spectrum observed by COMPTEL 
and INTEGRAL in Cygnus X-1. The main difference is 
caused by the energy injected into relativistic particles: the 
lower the hadronic content, the higher the power injected 
in relativistic particles. Neutrino production increases with 
a. 

A significant amount of flux in the range 10 MeV - 1 TeV 
is not expected because of absorption in the thermal photon 
field. We instead predict the existence of a bump at very 
high energies (-E > 1 TeV) . This high-energy emission may 
be detectable in the future from different sources. Since it 
is of hadronic origin, detections or upper limits can be used 
to place constraints on the number of relativistic protons 
in the corona. Orbital mod ulation can be impor tant owing 
to the variable absorption ()Romero et al .11201 Dal) . 

For sources where it is difficult to detect an electromag- 
netic flare, the neutrino production can nevertheless yield 
detectable events. Our results show that there may be in- 
stances in which a neutrino flare could be detected, but the 
ga mma-ray counterpart is not. 

iPetropoulou fc MastichiadisI (|2011l . I2012D studied the 
temporal signatures of electron and proton injection in a 
compact, magnetized source. Both photon quenching and 
synchrotron emission are responsible for the dynamical non- 
linear behavior of the system. Although these results can- 
not be directly extrapolated to our corona model, since 
interactions and background thermal photon fields cannot 
be ignored, future research may identify evidence of some 
kind of dynamical cycle in sources such as Cygnus X-1. 

6. Conclusions 

We have illustrated that a consistent treatment of the non- 
thermal emission from a magnetized corona can be imple- 
mented by solving the set of coupled differential equations 
for all particle species. 

Our application to Cygnus X-1 provided both a good 
fit to the observational data and interesting predictions for 
very high energy and neutrino instruments. 

In the future, we will explore flare episodes taking place 
in low-mass X-ray binaries. These are attractive objects 
for the application of our model, because in these systems 
the companion star has a weak radiation field. Since the 
absorption will be negligible, both the electromagnetic and 
the neutrino fiares may be detectable. 
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